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Abstract 

We solve the three-dimensional time-dependent Schrodinger equation for a few-cycle circularly 
polarized femtosecond laser pulse interacting with an oriented target exemplified by an Argon 
atom, initially in a Sp^, or 3pj^ state. The photoelectron momentum distributions show distinct 
signatures of the orbital structure of the initial state as well as the carrier-envelope phase of 
the applied pulse. Our ab initio results are compared with results obtained using the length-gauge 
strong-field approximation, which allows for a clear interpretation of the results in terms of classical 
physics. Furthermore, we show that ionization by a circularly polarized pulse completely maps out 
the angular nodal structure of the initial state, thus providing a potential tool for studying orbital 
symmetry in individual systems or during chemical reactions. 



PACS numbers: 32.80.Rm 



I. INTRODUCTION 



The recent decades have shown substantial progress in strong-field physics with pulsed 
lasers. Pulses with field strengths equivalent to the field of the Coulomb interaction in 
ground state atoms and durations of only a few femtoseconds are now available for a large 
range of wavelengths [H E] . This development has facilitated the opening of exciting new 
research areas such as attoscience [H |2] which, in turn has set entirely new standards for 
the interrogation of atomic and molecular dynamics. Strong-field ionization is a process of 
special relevance and importance within strong-field physics, since ionization triggers other 
strong-field phenomena, such as high-harmonic generation, the key process for generating 
coherent attosecond pulses in the XUV regime. In strong-field studies there is also an inter- 
est in creating unique quantum targets by molecular alignment and orientation techniques 
[3]. Spatial control of molecules is of importance within attoscience, femtochemistry and 
molecular reactivity among others. Recent experiments not only succeeded in alignment of 
molecules but also in 3D orientation, thus giving complete control over the spatial orientation 

The development described above motivates a detailed investigation of strong-field ion- 
ization of oriented targets by strong laser pulses with special emphasis on orbital structure, 
e.g., effects of orbital symmetry. The bulk of experimental and theoretical works within 
strong-field physics have dealt with linearly polarized laser pulses, and there is now a grow- 
ing interest in the study of ionization of atoms and molecules by elliptically polarized pulses 
[lOHTi] . The effects of angular orbital structure in the ionization by linearly polarized fields 
have previously been studied for direct electrons [151 IB] as well as for rescattered electrons 
[T7t [18] . When it comes to the probing of angular orbital structure through photoelectron 
angular distributions or momentum distributions, the use of a strong circularly polarized 
probe rather than a linearly polarized one has, however, two major advantages: (i) In the 
circularly polarized field, an electron born in the continuum is constantly driven away from 
the nucleus, due to the polarization. This dynamics minimizes rescattering effects as well as 
interference between wave packets launched at different instants of time during the driving 
pulse. Thus, the use of circularly polarized laser pulses entails a cleaner ionization signal 
with respect to orbital structure, (ii) The polarization plane is two-dimensional, which per- 
mits for a more transparent interrogation of angular orbital structure as we shall see in the 
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present work. In fact, it has recently been shown experimentally that strong- field ionization 
of 3D oriented C7H5N molecules by a circularly polarized field polarized in the nodal planes 
of the outermost orbitals provides a unique probe of the angular nodal structure [T^ . 

In this paper, we focus on the case where the laser polarization plane is perpendicular to 
the nodal plane. We calculate, by solving the time-dependent Schrodinger equation (TDSE), 
the photoelectron momentum distribution for ionization of an Argon atom initially prepared 
in a 3p^ or 3pj^ state modeling an oriented target with a single nodal plane. The and 3pj^ 
states, in particular, serve as models for investigating ionization of a molecular orbital with 
TT symmetry. In the strong near infrared laser field, the ionization is tunneling-like and the 
photoelectron is born in the continuum at a relatively large distance from the centre-of-mass. 
After the ionization, the circularly polarized field drives the electron away from the core, 
which in turn minimizes the importance of the detailed structure of the molecular potential 
at small distances. The calculations are compared with results obtained using the length 
gauge strong-field approximation (LG-SFA) [2ffl22], which facilitates interpretation using 
semiclassical theory [23]. Our results show distinct effects of the angular nodal structure 
of the initial orbital. For example the angular nodal structure of p^, p^^ and dxy orbitals is 
readily mapped out, showing that strong-field ionization by a circularly polarized laser pulse 
directly probes the nodal structure. 

This paper is organized as follows. In Sec. [Tl] we briefiy review the basic theory behind 



the computations. In Sec. Ill we discuss the results and in Sec. IV we conclude. Atomic 



units [h = Qq = rrie = 1] are used throughout, unless stated otherwise. 
II. THEORY 

We investigate an oriented target, modeled by an Ar atom initially prepared in a 3p^. or 
3pj^ state, interacting with a circularly polarized few-cycle pulse also in the xy plane. The 
electric field is defined as E = —dtA with the vector potential A{t), 



( cos [iijQt + 0) ^ 

sin (wot + 0) 

V " / 



(1) 



Here denotes the amplitude, cuq the carrier frequency, the carrier-envelope phase (CEP) 
and /(t) = sin^ (|^) the envelope, with the number of optical cycles. The laser field 
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couples primarily to the least bound electron, which in turn motivates the use of the single- 
active electron approximation (SAE). Hence, our analysis is restricted to the outermost 
electron, the remaining being described by an effective potential [21j. The photoelectron 
momentum distribution is given by 

^^-g-H(*?«l*(r^,r))^ (2) 

where ^E'^(r) is a continuum scattering wave function of asymptotic momentum q and \E'(r, T) 
is the wave packet at the end of the pulse, T. We calculate the wave packet by expand- 
ing in a spherical harmonic basis for the angular part and reduced radial wave functions 
fim{r) defined on a radial grid, ^(f,t) = Y!i=o Y!m=-iUim{r,t)/r)Yim{9,4>), and then solv- 
ing the TDSE in the velocity gauge using a split-step method. The details are described 
elsewhere |25l |26]. The scattering states entering in Eq. (2) are obtained by solving the 
time- independent Schrodinger equation for the effective potential. 

The LG-SFA [201121], which completely neglects the Coulomb interaction in the final state 
as well as all intermediate states, offers a simple formula for the momentum distribution 

8^P r'^ ^ 

-—^^l^^^v^rME-A^mm . (3) 

where is a Volkov wave function with asymptotic momentum q and \E'j is the ini- 
tial state. In the evaluation of Eq. (3) it is accurate to use the asymptotic form, 
\E'j(r, t) = Qm'"'^"^ exp(— Kr)y;m(^5 0) exp(z/pt), of the initial state, since ionization pri- 
marily occurs at large distances from the nucleus. Here Ip denotes the ionization potential, 
K = \/'il-p and V = ZjK with Z the charge of the residual ion. The spatial integration in Eq. 
(3) can be performed analytically, within this approximation, while the ID time integral 
is evaluated effectively using the saddle-point approximation [27j, which is accurate in the 
parameter regime considered in this paper [28| 129] . 

III. RESULTS 

A. Imprints of nodal planes in momentum distributions 

Figure 1 shows the TDSE results for the momentum distributions, in the xy plane of 
polarization, for the (a, b) 3p^ and (c, d) 3py initial states of Ar, probed by an 800 nm 
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3-cycle laser pulse with peak intensity / = 1.06 x lO^^W/cm^ and CEP (a, c) (f = — 7r/2 , 
(b, d) Lp = 0. The distributions are calculated using a 4096 points radial grid extending to 
'"max = 400, maximum angular momentum /max = 40, and a time-step of 0.005. Notice that 
the ionization potential for the Ar(3p) state is 15.76 eV. The low-energy part a/q'x- + Qy < 0.1 
has been removed for better graphical display. For the Sp^, state with = — vr/2. Fig. 1(a), 
and 3p^ state with = 0, Fig. 1(d), we observe a distribution with a single dominant peak, 
which, as we shall discuss in more detail below, is located near q ^ —A{T/2). In the two 
other cases we observe a splitting of this peak into two nearly symmetric peaks. As we now 
show these main features are explained using the LG-SFA. Figure 2 shows the momentum 
distributions obtained using the LG-SFA, in the plane of polarization, for the (a, b) 3p^ and 
(c, d) 3p^ initial states of Ar, at 800 nm light and peak intensity / = 1.06 x lO^^W/cm^, 
CEP (a, c) ip = — 7r/2, (b, d) ip = and three optical cycles = 3. 

The SFA results clearly reproduce the general structures of the TDSE distributions, 
except for an overall angular shift and a spiral like interference structure. The former overall 
angular shift is due partly to the Columbic interaction, partly to the fact that the laser field 
breaks the rotational invariance perpendicular to the polarization plane [l^. This feature 
has been observed experimentally in ionization of Helium by a few-cycle circularly polarized 
laser pulse [llj. The spiral structure is connected to the interaction between outgoing low- 
energy electrons and the long-range Coulomb potential, equivalent to the radial fans seen in 
the linear case [HOl - IM] . 

It is well-known, in the case of a spherically symmetric initial state, that the LG-SFA 
momentum distribution reflects the shape of the vector potential [28] . This is easily explained 
by semiclassical theory [23]: Due to the exponential dependence of the rate on the field 
strength, ionization predominantly occurs at the times t when the magnitude of the electric 
field is close to the maximum \E{T/2)\. Assuming that the (classical) electron is born in the 
continuum with zero initial velocity and that it subsequently moves under the influence of the 
electric field only, we obtain ggnai = —A{ti) with ti the instant of ionization (ggnai = —A{T/2) 
dominates). However, in the case of a non-spherical initial state, this simple argument has 
to be extended to take into account the shape of the initial state, in particular the angular 
nodal structure. For example, it may happen that the field is in a nodal plane of the initial 
state at t = T/2 and hence ionization will be suppressed. In other words, it is in this more 
general case, the shape of the vector potential, electric field and initial state that determines 
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FIG. 1: Momentum distributions from the single-active-electron TDSE calculation, for the (a, b) 
3p^ and (c, d) Sp^^ states of Ar, in the xy polarization plane, for the following choice of laser 
parameters: Angular frequency u = 0.057, corresponding to 800 nm, peak intensity / = 1.06 x 
lO^^W/cm^, carrier envelope phase (a, c) </? = —7^/2, (b, d) ip = and number of optical cycles 
N = 3. 
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FIG. 2: SFA momentum distributions, for the (a, b) 3p^ and (c, d) 3py states of Ar, in the xy 
polarization plane, for the following choice of laser parameters: Angular frequency u = 0.057, 
corresponding to 800 nm, peak intensity / = 1.06 x lO^^W/cm^, CEP (a, c) ip = -7r/2, (b, d) 
(/? = and number of optical cycles = 3. The white curves show a parametric plot of —A{t) 
from Eq. (1). If no nodal planes are along the polarization at time t = T/2, the dominating peak 
in the distribution follows the simple formula gJinai — —A{T/2) as seen in (a) and (d). If on the 
other hand the field peaks in the direction of the angular node (b), (c), the dominating peak in 
the distribution splits, simply because there is no charge density in that direction. 

7 




-^3 



b 




mi 







x(a.u.) 



FIG. 3: Electron density in the Sp^, state of Ar together with a parametric plot of the electric field, 
corresponding to the vector potential shown in Fig. 2, scaled by a factor. In (a), (p = —it/2, the 
electric field at the peak of the pulse t = T/2 points in a direction where the orbital density is 
largest. In (b), 99 = 0, the electric field at the peak of the pulse runs through the orbital node (no 
electron density). 



the overall structure of the momentum distribution. The results in Figs. 1 and 2 clearly 
illustrate this combined effect. Let us start by analyzing the results for the 3p^. state and let 
us concentrate on the LG-SFA results to make the discussion as simple as possible. Figure 
2(a) shows the SFA momentum distribution in the plane of polarization for = —tt/2 and 
3p^ initial state. We observe a distribution with a single peak located near —A{T/2), in 
good agreement with the semiclassical theory. In the case of (p = —tt/2, E{T/2) = Ex{T/2)x 
points in the direction of maximum electron density (see Fig. 3(a)), while both E and 
are symmetric with respect to {x,y) — )■ {x,—y). This leads to a single peak situated at 
q ~ —A(T/2) oc —y and symmetric with respect to {x,y) — )■ {—x,y) in good agreement 
with the calculated distribution. Figure 2(b) shows the SFA momentum distribution, in the 
plane of polarization, for = and a 3p^ initial state. We observe two peaks, symmetric 
with respect to {qx,<iy) — > {(lx,—%), in contrast to the simplest single peak semiclassical 
prediction. In the case of = 0, E{T/2) = Ey(T/2)y lies in the nodal plane of 3p^ (see 
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Fig. 3(b)), i.e., there is no electron density to ionize ect t — T/2 and hence we do not have 
ionization at field maximum. When \t — T/2\ increases from zero the electric field begins 
to point in directions with non-zero electronic density. Thus, the single peak from before 
splits into two symmetric peaks, due to the fact that E and Sp^, are symmetric with respect 
to {x,y) — )■ {—x,y), again in good agreement with the obtained distribution. The precise 
location of the peaks is determined by a competition between the rate arising from the 
electric field and the density profile of the initial state. The results shown in Fig. 2(c) and 
Fig. 2(d) are explained using similar reasoning. Notice, that the signature of the nodal 
plane in the momentum distribution is advanced by 7r/2 compared to the orbital angular 
nodal structure, reflecting the n/2 phase between the electric flcld and the vector potential. 

The above discussion of Figs. 1, 2 and 3 clearly shows that the nodal structure of the 
initial orbital is mapped uniquely to the momentum distribution when the laser polarization 
plane is perpendicular to the nodal plane. This effect is even more visible in the case of a 
shghtly longer femtosecond pulse, where the asymmetry caused by the CEP of the pulse is 
smaller. Figure 4 shows the SFA reduced momentum distribution in the plane of polarization 
d'^P/dqxdqy = J {d^P /dqxdqydqz)dqz for H(ls), H(2p2,), H(2pj^) and H(3da;j^) states obtained 
using the same parameters as in Figs. 1 and 2 except that the number of optical cycles was 
increased io N — 2Q and / = 1.0 x lO^^W/cm^. The use of the LG-SFA is for computational 
convenience and is justified by the results presented on Ar. The binding energy of the 
H(2p^), H(2p^) and H(3dj;j/) have been artificial modified to —0.5, by modifying k — ^J2T^^ 
in order to compare with the H(ls) result. Figure 4 clearly shows that the angular nodal 
structure is mapped uniquely to the momentum distribution. This is particular clear when 
compared with Fig. 5, which shows the asymptotic densities for the Is, 2p^, 2p,y and ?)dixy 
states. The electric field is almost constant during the dominant cycles of the field, i.e., 
there is no competition between the rate arising from the electric field, and the density 
profile. In the case of H(ls) we observe a completely symmetric distribution with well- 
resolved above-threshold ionization peaks, while ionization is strongly suppressed along the 
X-axis (|/-axis) in the case of H(2p^) (H(2pj^)) and along the x-axis and ?/-axis in the case 
of H(3da;j^), reflecting the angular nodal structures in Fig. 5. Ionization is suppressed along 
the X-axis (y-axis) in the case of H(2p2,) (H(2pj^)), due to the yz [xz] nodal plane which is 
shifted 90° by the vector potential. There are two angular nodal planes for H(3dj;2^), namely 
the xz and yz planes, suppressing ionization along the y-axis and x-axis. 
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FIG. 4: SFA momentum distributions d'^P/dqxdqy = f{d^P/dqxdqydqz)dqz for (a) Is, (b) 2p^, 
(c) 2py and (d) Sd^^j^ obtained using the following laser parameters: Angular frequency oj = 0.057, 
corresponding to 800 nm, peak intensity / = 1.0 x lO^^W/cm^, carrier envelope phase = 
and number of optical cycles N = 20. The binding energy of the 2p^, 2py and Sd^y states are 
modified to the ground state energy of H. The white dashed lines show regions where ionization is 
suppressed. 
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B. Carrier-envelope phase effects in circularly polarized fields 



In Figs. 1 and 2, we note that changing the CEP from = — 7r/2 to ip = does not 
correspond to a counterclockwise rotation of angle 7r/2, in contrary to what one might expect: 
CEP effects for atoms and linear molecules aligned along the propagation axis of the applied 
circularly polarized field were characterized in terms of simple rotations of the total system 
around the propagation axis [35J. This characterization, however, only holds when the initial 
state is invariant with respect to rotations around the propagation direction or when the 
system initially is described by a uniform incoherent mixture of magnetic substates. The 
latter point was not discussed exphcitly in [55] and we use the opportunity to do so now. 
Note that CEP effects have been described generally in a Floquet-like approach [36l |37] . 

Consider an atom or a linear molecule, described within the SAE, interacting with a 
laser pulse described by the vector potential A{t) defined in Eq. (1). We assume that the 
field-free Hamiltonian Hq is invariant under rotations around the 2;-axis. Furthermore, we 
assume that the probability for a specific magnetic quantum number Pm is uniform. The 
ensemble average of d^P/dq^dgydgz is defined as 



Tr(p(0;t)P, 



<?/' 



(4) 



dq^dqydq, 

where t is any time after the pulse, p(0; t) is the density matrix of the total system, 
the CEP and P^- = |\E'Z)(\EfZ| projects on the exact scattering states |\E'^) with asymptotic 
momentum q. The density matrix is given by 



nJM 



(5) 



M=-J 



Here Pm = 2J+1 |^njM(0;i)) = U{(f);t,0)\nJM) , where U is the time-evolution oper- 
ator for the total system and \nJM) the field-free initial state with J denoting the total 
angular momentum of the system, M the corresponding magnetic quantum number and n 
the remaining quantum numbers. We note, in passing, that the Hamiltonian and hence the 
time-evolution operator have a parametric dependence on (p through the interaction with 
the external field. We evaluate the ensemble average Eq. (4) using the position-eigenstate 
basis and thereby obtain 



Q3p 



dq^dqydq, 



2J + 



1 ^ 



M=-J 



'nJM\ 



(6) 



11 



Using the relation j39] 



^^^m'm("' 7)^m'*m(«, /3, 7) 



-Sm'm" 



2J+1 

where D^^,j^j{a, /3, 7) are the Wigner rotation functions, (a, /3, 7) the Euler angles and 
sin{P) d (3 dad'y, Eq. (6) can also be expressed as 



(7) 



Q3p 



dqxdqydq^ 



■in" 



dVt 



df{mZy^!nJAp{^A\t) 



(8) 



Here \E'„jAf'(fi, 0; t) is a solution to the TDSE corresponding to the rotated initial state 
D{VL)\nJM'), with D the rotation operator and = (a,/3,7). However '$nJM'{^,4'',t) = 
exp(-zJ^0)^„jM'(^',0 = 0;t) with fi' = (a -0,/3,7) [33]. Thus 



Q3p 



dq^dqydq^ 



in" 



dQ 



(9) 



The rotation of a scattering wave function with asymptotic momentum q can be ac- 
complished just by rotating the asymptotic momentum, i.e., exp{iJz(f))'^'^ = \E'^,, where 
q' = Iizi—<P)q, with the 3x3 orthogonal matrix which generates counterclockwise rota- 
tions around the 2;-axis [39 1 . This means that 



Q3p 



dq^dqydq^ 



1 

8^ 



2n c2-K 



^0 



sin(/3)(i/3(ia(i7 



y dr{^--.X^r.JM'{^A = ^\t) 



(10) 



Finally this expression can be rewritten as 



Q3p 



dqxdqydqz 



1 



in" 



2-K /.27r 



sin(/3)(i/3(ia(i7 



Q3p 



-(g^,0 = O) 



(11) 



dq^dqydq^ 

due to the uniformity of the distribution function over the orientations, G{Q) = l/{8n'^). 
Thus a change in the CEP from = to = 0' corresponds to a rotation of the system 
around the 2;-axis by 0'. This is illustrated in Fig. 6, which shows the ensemble average 
of the exact momentum distribution, in the xy polarization plane, for Ar atoms initially in 
an incoherent mixture of 3p^. or 3py states, for two different values of the CEP = —n/2 
and = 0. Notice, that the validity of Eq. (11) does not require the SAE or the dipole 
approximation as can easily be checked by going through the steps leading to Eq. (11) 
maintaining full retardation and accounting for all electrons. 
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The pure Sp^. and Sp^^ states of Ar, of course, do not fulfil the requirements for the 
validity of Eq. (11) and hence the CEP effects observed in this paper are not in general 
simple rotations. It is, nevertheless, still possible to express momentum distributions for Sp^, 
in terms of momentum distributions for Sp^. Symmetry considerations immediately leads 
to the following formula 

^3 p 03 p 

^-^ + vr/2) = R,(7r/2)— ^(0), (12) 



dq^dqydq^ dq^dqydq^ 

where (f) is the CEP. This geometric correspondence is illustrated in Figs. 1 and 2, where 
we note that d^Ps^ /{dq^dqydq,){^ = 0) = R,{n/2)d^P3pJ{dq^dqydq,){(f) = -n/2). 



IV. CONCLUSION 

We have investigated strong-field ionization of Ar(3p^.) and Ar(3py) in the presence of a 
strong circularly polarized laser pulse, by calculating the photoelectron momentum distri- 
bution using the TDSE. These systems model the ionization from an oriented target with a 
single angular nodal plane interacting with a circularly polarized laser pulse, a scenario that 
could be realized in an oriented molecular target. Our results are compared with results 
obtained using the LG-SFA and, up to a relatively small rotation of the distribution that is 
well-understood as a combined field and potential effect [H], good agreement is observed, 
emphasizing the LG-SFA as a valuable tool for studying ionization by circularly polarized 
pulses. Both TDSE and LG-SFA results show distinct signatures of the initial states as well 
as the temporal shape of the applied pulse. The findings indicate that strong- field ionization 
by circularly polarized pulses could serve as a probe for revealing nodal structure and also 
changes that may occur during chemical reactions. 
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FIG. 5: The asymptotic density |Eime"'"^mP (a) H(ls), (b) H(2p^), (c) H(2pj^) and (d) 
il{3dxy The white dashed hnes show the angular nodal planes. 



16 




q^(a.u.) q^Cau.) 

FIG. 6: Ensemble average of the SAE, TDSE momentum distributions in the xy polarization plane 
for Argon initially in a uniform incoherent mixture of Sp^, and 3pj^ states. The laser parameters 
are: Angular frequency lo = 0.057, corresponding to 800 nm, peak intensity / = 1.06 x lO^^W/cm^, 
carrier envelope phase (a) (p = —it/2, (b) (f) = and number of optical cycles = 3. 
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